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ABSTRACT 

We investigate the use of simple colour cuts applied to the Sloan Digital Sky Survey (SDSS) 
, optical imaging to perform photometric selections of emission line galaxies out to z < 1. Our 

Qh' selection is aimed at discerning three separate redshift ranges: 0.2 ^ z < 0.4, 0.4 ^ z < 0.6 

and 0.6 ^ z < 1.0. We calibrate the selections using data taken by the COMBO-17 survey 
^ , in a single field (Sll), which is covered by the SDSS imaging. We thus perform colour-cuts 

^ ' using the SDSS g, r and i bands and obtain mean photometric redshifts of z/oto — 0.32 ± 0.08, 

I Zmid — 0.44±0.12 and z/ii = 0.65±0. 21. We further calibrate our high redshift selection us- 

ing spectroscopic observations with the AAOmega spectrograph on the 4m Anglo- Australian 
Telescope (AAT), observing « 50 — 200 galaxy candidates in 4 separate fields. With just 1- 
hour of integration time and with seeing of « 1.6", we successfully determined redshifts for 
« 65% of the targeted candidates. We compare our spectroscopic redshifts to the photometric 
I redshifts from the COMBO-17 survey and find reasonable agreement between the two. We 

(•~^ , calculate the angular correlation functions of these samples and find correlation lengths of 

ro = 2.64i°:°^/i-iAfpc, tq = 3.62 ± 0.06h-^Mpc and vq = 5.88 ± 0.12h-^Mpc for the 
low, mid and high redshift samples respectively. Comparing these results with predicted dark 
matter clustering, we estimate the bias parameter for each sample to be 6 = 0.70, b = 0.92 
' and b = 1.46. We calculate the 2-point redshift-space correlation function at z « 0.6 and find 

I a clustering amplitude of Sq ~ 6.4 ± Q.8h^^ AIpc. Finally, we use our photometric sample 

to search for the Integrated Sachs-Wolfe signal in the WMAP Syr data. We cross-correlate 
our three redshift samples with the WMAP W, V, Q and K bands and find an overall trend 
for a positive signal similar to that expected from models. However, the signal in each is rel- 
■ atively weak with the results in the WMAP W-band of WTg{< 100') = 0.25 ± 0.27 fiK, 

I 0.17 ± 0.20//A' and 0.17 ± O.lQiiK for the low, mid and high redshift samples respectively. 

Combining all three galaxy samples we find a signal of WTg{< 100') = 0.20 ± 0.12pK in 
the WMAP W-band, a significance of 1.7a. However, in testing for systematics where the 
WMAP data is rotated with respect to the ELG sample, we found similar results at several 
different rotation angles, implying the apparent signal may be produced by systematic effects. 

Key words: cosmology: photometric redshifts, clustering Integrated Sachs-Wolfe Effect. 
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1 INTRODUCTION 

Imaging surveys are currently in the process of mapping out 
a vast region of the Universe over a large range of the elec- 
tromagnetic spectrum. The pace- setter in recent years is the 
Sloan Digital Sky Su rvey (SDSS, lYork et alj|2000l). which now 
provides (as of DR6, lAdelman-McCarthv et al. I l2008h photomet- 
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ric data for approximately 230 million distinct sources over an 
area of 8240 deg^ . Curr ent and future wide and deep field 
surveys such as SWIRE (McMahon et al. 2001 ; Irwin & Lewia 
2001 ; Rowan-Robinson et al. 2008), UKIDSS ( Lawren ce et all 
nOOi) . CFHTLS (Caban ac et al.l l2007l: iMcCracken et alJ l2008h . 
Pan-STARRS tKaiseru200l) . LSST l lSweenev et all 1 20091) . VST 
ATLAS, the VISTA Hemisphere Survey and the SDSS itself will 
continue to add to the mapping of the Universe around us present- 
ing increasing amounts of data at a variety of wavelengths. Given 
this enormous effort in the collection of photometric data, and the 
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expense of subsequent spectroscopic surveys, the filtering of galax- 
ies by type and redshift via their photometric properties is a vital 
and powerful tool for the effective use of the large quantities of 
photometric data available to us. Selecting distinct galaxy popu- 
lations in this way offers a relatively cheap route to large galaxy 
and QSO surveys, either through using photometric redshifts or by 
using broader photometric constraints to select specific galaxy pop- 
ulations for subsequent spectroscopic surveys. 

A key example of the success of this process over re- 
cent years has been the photometric selec tion of Luminous Red 
Galaxies (LRGs) at redshifts up to z~0.8 (Eisenstein et al. 2001; 
IPadmanabhan et alj20053 ; lBlake et alTlOOV. : Collister et al. 2007) . 
LRGs offer a simple route to photometric selection of different 
redshift samples due to the 4000A break feature, which passes 
through the optical wavelength bands from redshifts of z = 
out to z = 1. Such selections have been used in a range of key 
cosmological applications and perhaps most significant amongst 
these is the detection of the baryon acoustic oscillation (BAO) sig- 
nature in analyses of large scale structure. This was measured by 
[Eisenstein et al. (2005) using spectroscopic redshifts of > 40, 000 
LRGs selected from the SDSS using photometric constraints based 
on the progression of the 4000A break. The measurement of 
the BAO signal should also be possible using suitably accurate 
photometric redshifts as the B AO galaxy clustering features are 
on sc a les of order 100h~^M pc jEisenstein et alj 200 ll: Cole et al 
I2OO5I: lEisenstein et alj |2005|) . In this vein ISawangwit & Shanks 
1 I2OO9 I) used some 500 LRGs with spectroscopic redshifts obtained 
using AAOmega on the AAT to calibrate an r — i ; i — z colour-cut 
in terms of redshift and obtained a redshift distribution in the range 
0.55 < z < 1.0, with z — 0.68±0.1. They go on to use a full pho- 
tometric sample based on this selection, combined with two others 
with z « 0.35 and z ^ 0.55, to investigate the BAO signal. Us- 
ing these photometric selections to measure the BAO signature in 
the angular correlation function, they sh ow an agreement with the 
spectroscopic measurements reported by Eisenstein et 'Zl( l2005[) . 

A further example of the use of photometrically selected 
LRGs is the study of the Integrated Sachs Wolfe effect (ISW), 
in which CMB photons are subjected to shifts in energy as they 
pass through the gravitational pot ential wells of galaxie s and clus- 
ters in an accelerating Universe dSachs & Wolfein"967h . This has 
been studied by a number of authors using the WMAP 1st and 
3rd year data releases in combination with photometrically selected 
galaxy populations. Asi de from the use of simple magnitude lim- 
ited galaxy samples (e.g.lposalba et al .l2003l : lFosalba & Gaztanag3 



l2004l ; lRassatetal.ll2007|)^ this work is again domin ated by the use 
of LRGs. For example. IPadmanabhan et al] ( l2005bl) use photomet- 
rically selected LRGs from the SDSS (covering a redshift range 
of 0.2 < z < 0.6) to make a 2.5g detection o f the ISW effect in 
WMAP 1 St year data, whilst lCabre et alj ( l2006h combine a 2 ~ 0.5 
SDSS LRG sample with the WMAP 3 rd year data to claim a de- 
tected ISW signal of 3a. More recentlv rGiannantonio et al] ( l2008h 
performed a combined analysis using a number of galaxy a nd QSO 
catalogues, incorporating 2MASS data 1 Jarrett et al. 200Q) . SDSS 



DR6 galaxi es, SDSS MegaZ L RGs iCoUister et al j|2007h. NVSS 
radio data jCondon et aklligQSh . HEAO X-ray data ('Boldt"l987') 
and the SDSS DR6 QSO catalogue feichards et al. 2009). By com- 
bining the results from these datasets, iGiannantonio et alj ( |2008|) 
report a 4.5a detection of the ISW effect and go on to use this to 
test a number of cosmological models. In particular, their ISW re- 
sult places constraints on the mass density of n„i = 0.26lo;o?, 
assuming a flat ACDM cosmology. 

As an alternative to the recent dominance of LRGs for large 



redshift surveys, we now look at the photometric selection and clus- 
tering of emission line galaxies (ELGs). The key advantage of using 
ELGs in this application is the ability to identify galaxies through 
spectroscopic observations with relatively short exposure times, 
due to the emissio n lines with which they are most easily identified. 
Blak e et al] ( 1200^ have successfully used this to their advantage to 
undertake the WiggleZ spectroscopic survey of ELGs with the aim 
of measuring the BAO signal at z«0.7. They select ELGs in the red- 
shift range 0.5 < z < 1 using a combination of GALEX FUV and 
NUV imaging with g, r and i band optical imagin g from SDSS an d 
the second Red Cluster Sequence (RCS2) project ( lYee et alj2007h . 
In this case the GALEX data allows them to select z > 0.5 galax- 
ies using the Lyman Break technique, whilst the SDSS optical data 
allows them to limit their sample to just blue emission line objects. 
They then perform spectroscopic observations using the AAOmega 
instrument at the 4m AAT, with exposure times of < Ihr required 
for successful identification. 

In this paper we attempt to refine a number of photomet- 
ric selections of ELGs in different redshift ranges, using optical 
data from the SDSS. In making greater use of the ELG popula- 
tion in studies of large scale structure, we may maximise the use 
of the available data from such large scale surveys as the SDSS 
and the upcoming VST ATLAS and Pan-STARRS surveys. We use 
COMBO- 17 photometric redshift data in combination with SDSS 
data to perform a calibration of our photometric selections (sec- 
tion |2j. We go on to outline and review our observation program 
at the AAT, which was aimed at providing a spectroscopic cata- 
logue with which to further calibrate the photometric selections. In 
section [3] we then evaluate the clustering properties of the galaxy 
populations contained in our photometric selections using SDSS 
data. We then use the full samples of > 600,000 galaxies selected 
from the SDSS to perform a search for the ISW effect in WMAP 
5yr data (section |4). We present our summary and conclusions in 
section |5] Throughout this paper we assume a ACDM cosmology 
with = 0.3, = 0.7 and Ho = WOhkms'^Mpc. 



2 PHOTOMETRIC SELECTION 
2.1 Data & Selection 

Using SDSS imaging data, our aim is to develop a set of photomet- 
ric selection criteria using the SDSS filter bands tFuku g ita 19961) 
alone to isolate emission line galaxies in 3 separate redshift ranges 
of approximately z < 0.4, 0.4 < z < 0.6 and z > 0.6. At these 
redshifts, the 4000A break is a key feature in the observed optical 
spectra of both red and blue galaxies as it moves through the g and 
r SDSS filters with increasing redshift. In ELGs however, the break 
is somewhat weaker than in the spectra of LRGs, whilst the contin- 
uum at wavelengths greater than 4000A remains lower compared to 
the LRG spectrum due to the dominance of young blue stars in the 
ELGs. These contrasts in the spectra of LRGs and ELGs inherently 
allow us to separate the two in colour space, whilst simultaneously 
facilitating photometric selections of galaxies at different redshifts. 

With this in mind we have used the Bruzual and Chariot stellar 
population synthesis code (Bruzual & Chariot 2003) to model the 
evolution of a typical emission line galaxy in the gri (AB) colour 
plane. We used a Salpeter IMF with a galaxy formed at z = 6.2 
(i.e. with an age of 12.6Gyr at 2 = 0) and a r = 9Gyr exponen- 
tial SFR. The resultant gri colour evolution track from z — 1.2 
to 2 = is shown in figure [T] (dashed black line). Here we see a 
clear evolution in the gri colour space around which we may build 
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Figure 1. Galaxies in the COMBO-17 Sll field plotted in the gri colour 
plane using SDSS magnitudes. The cyan diamonds, green tiiangles and red 
squares show galaxies with photometric redshifts in the ranges 0.2 < z < 
0.4, 0.4 < z < 0.6 and z > 0.6 respectively. The labelled black tracks 
show the evolution of an emission line and an elliptical galaxy based on a 
Salpeter model, from z = 1.2 to z = 0. Our photometric selections are 
mai'ked by the dotted blue, solid green and dash-double-dotted boxes for 
our low, mid and high redshift bins. 



a selection regime for identifying candidates in our desired redshift 
ranges. We also plot a track (dot-dash line) for an elliptical galaxy 
using a T = IGyr exponential SFR and a redshift of formation of 
z=10 (with Solar metallicity and Salpeter IMF). 

We calibrate our selections using the photometric redshift data 
publi shed by the COMBO-17 team CWolf et alJi 20d3l. ISimon et al.l 
|2008|) . The data we use is from the COMBO-17 Sll field, which 
covers an area of 0.5° x 0.5° centred at 1 Ih 42m 58s, -01 42' 50" 
(J2000) and the entirety of which is covered by SDSS imaging. 
These data provides accurate (dz/{l + z) — 0.02) photometric 
redshifts for a total of 7248 galaxies based on broad and narrow 
band imaging using 17 different optical and near-infrared filters. 
We match the positions of COMBO-17 galaxies to the equivalent 
objects in the SDSS data, thus combining the SDSS ugriz magni- 
tudes with the COMBO-17 photometric redshifts. The Sll objects 
are shown in the (SDSS) g — r, r — i colour plane in figure[T] For 
the purposes of clarity, we only plot those galaxies classed as blue 
spirals by the COMBO-17 team in this figure, however the location 
of red-sequence galaxies is indicated by the LRG evolution track 
(dot-dash line). The COMBO-17 galaxies have been split into three 
populations for the purposes of this plot based on their assigned 
photometric redshift from COMBO-17: 0.2 < z < 0.4 (blue dia- 
monds), 0.4 < z < 0.6 (green triangles) and 0.6 < z < 1.0 (red 
squares). 

Based on the distribution of the photometric redshifts and the 
ELG evolution track presented in the above plot, there is a clear 
progression in the gri colour plane based on ELG redshift. Further 
to this, areas of the plot can be isolated that should minimize the 
number of red-sequence galaxies, whilst maximizing the numbers 
of either z < 0.4 or 2 > 0.6 ELGs. The medium redshift range 
does however present significant problems. The ELG evolution 
track appears to pass through a region populated by both lower and 



Figure 2. Redshift distributions of our three photometric selections based 
on photometric redshifts from the COMBO-17 data in the Sll field. The 
three samples give mean redshifts of 2;o = 0.29 it 0.05, Zmid = 0.44 it 
0.08 and z,ji = 0.65 it 0.21. 



higher redshift ELGs as well as low-redshift red-sequence galaxies 
in the 0.4 < z < 0.6 range. 

From the above observations we construct three sets of colour 
cuts to preferentially select three redshift ranges. These are shown 
in figure[T]by the dotted blue box (low redshift cut), solid green box 
(medium redshift cut) and dash-double-dot red box (high redshift 
cut). As discussed above, the mid redshift range is significantly 
exposed to contamination from both ELGs at unwanted redshifts 
and red-sequence galaxies. To minimize the numbers of these we 
have therefore added colour cuts to this sample based on the r ~ i, 
i ~ z and u — g colours of the selected galaxies. These additional 
cuts have also been calibrated using the COMBO-17 photometric 
redshifts. The details of our three selections, including the addi- 
tional mid-redshift colour cuts, are given explicitly in table[T] These 
cuts have been tailored to produce sky densities of candidates of 
~ IQQdeg^^ for each of the three redshift ranges in order to pro- 
vide candidate numbers suitable for wide field spectroscopic sur- 
veys performed with instruments such as the 2dF/AA0mega spec- 
trograph. 

The photometric redshift distribution for our three samples are 
shown in figure |2] This plot includes all the selected galaxies from 
the Sll field, including those identified as being part of the red- 
sequence (these making up ~ 4% of the total selected across all 
three selections). The three selections are characterised by mean 
redshifts of ziow ~ 0.29, Zmid = 0.44 and zm — 0.65, with 
standard deviations of aiow ~ 0.05, aiow = 0.08 and ahi ~ 0.21. 

For the purposes of this paper we now use our three selections 
to create three datasets from the SDSS galaxy catalogue. We ap- 
ply the selections to the SDSS DR6 data release taking our data 
from the PhotObjAll table. Aside from the colour-magnitude crite- 
ria given in table[T] we reject objects which do not meet the follow- 
ing criteria: 

• TYPE = 3 (i.e. classed as a galaxy); 

• NCHILD = 0; 

• Flagged as BINNED 1, BINNED2 and BINNED3; 

• 90° < R.A. < 270°. 
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Table 1. Selection criteria chosen to identify galaxies in our three redshift ranges: 0.2 < z < 0.4, 0.4 < z < 0.6 and 0.6 < z < 1.0 using SDSS ugriz AB 
magnitudes. These are illustrated in the ugr colour plane in figure[T] 



Low redshift Mid-redshift High redshift 

19.0 < i < 20.0 19.0 < i < 20.2 19.5 < i < 20.5 

r - i < 0.3 g - r > 1.2(r - j) + 0.1 r - j > 0.5 

g - r < 1.2(r - i) + 0.9 g - r < 1.2(r - j) + 0.9 3 - r < 1.2(t- - i) + 0.06 

9 - r > -1.2(r - j) + 0.75 g - r > -1.2(r - i) + 1.65 3 - r > 1.2(r - j) - 0.6 

g - r < -1.2{r - i) + 1.3 g - r < -1.2{r - j) + 1-95 g ~ r < -1.2{r - i) + 2.2 

r - i > -{i - z) + 0.5 
-2.0 <u-g < 1.0 
i - z> 0.55 



We also limit our selection to the main SDSS region, rejecting 
stripes 40-43. Stripe 26 is also rejected as this appears to show 
some contamination and spurious density fluctuations. The photo- 
metric selections are performed using the SDSS model magnitudes 
with the appropriate extinction values subtracted. The total num- 
bers of candidates given by each selection are 892,528, 620,020 
and 734,566 for the low-, mid- and high redshift selections respec- 
tively. These numbers give sky densities of W3deg~^, 71.9deg~^ 
andSS.ldeg"^. 



2.2 Observations 

An important element of this work is the calibration of the photo- 
metric selection samples with spectroscopic observations to con- 
firm the achievable redshift distribution of our selections. To 
this end, we have performed spectroscopic observations of our 
z g:i 0.7 sample usi n g the AAOmega spectrograph at the AAT 
dSaunderl et all |2004 rSharp et al.l I2OOT ). AAOmega is a double 
beam spectrograph fed by 2.1" diameter fibres, which allows the 
simultaneous observations of up to ~ 360 objects in a circular field 
of view of diameter 2°. 

Observations were taken on the AAOmega instrument at the 
Anglo- Australian Observatory (AAO) on the nights of 4th and 6th 
March 2006. The spectrograph was configured using the 5700A 
dichroic, with the 580V grating mounted in the blue-arm and the 
385R grating in the red-arm. The 580V grating gives a wavelength 
coverage of 370nm to 580nm, with a pixel size of O.lnm/pixel 
and the 385R a coverage of 560nm to 880nm, with a pixel size 
of 0. 16nm/pixel. Both provide a resolution of 1300. In total, 
AAOmega offers 400 fibres per observation, however a signifi- 
cant number of these were at times used for other projects (e.g. 
IRoss et alj200i) . were locked to guide-stars, sky-targets or simply 
malfunctioning and so our target numbers range from ~50-230 per 
observation. We targeted four 2° fields with multiple exposures of 
1200s each. The observations are summarised in table |2] The four 
observed fields are targeted towards the Cosmic Evolution Survey 
(COSMOS, Scoville et al. 2007) field, the d05 and e04 fields from 
the 2dF-SDSS LRG and QSO (2SLA0, ICroom et all2004h survey 
and the SI 1 field from COMBO-17. 

Target objects were selected using our selection criteria ap- 
plied to the SDSS data available for each of the fields. Over the 
course of the 5 nights, seeing ranged from ~ 1.5" — 3.0", with a 
mean of ~ 2.0". All observations were flat fielded, arc calibrated 
and combined using the AAO's 2dFDR tool. Approximately 20% 
of fibres were affected by an early instrumentation problem known 
as fringing, which led to an almost sinusoidal signal in the output. 
In the SI 1 field this affected 27 of the fibres targeted on ELG can- 
didates. A further problem, is encountered due to the strong sky 



Table 2. Co-ordinates of the four fields targeted with the number of gri 
selected ELG candidates in each 2-degree field. 



Field 


COSMOS 


d05 


e04 


Sll 


Date 


06/03/06 


04/03/06 


04/03/06 


06/03/06 


R.A. 


150.118° 


200.399° 


221.899° 


175.741° 


Declination 


+2.2052° 


-0.2124° 


-0.2141° 


-1.7159° 


Exposure time 


3 X 1200s 


3 X 1200s 


4 X 1200s 


3 X 1200s 


Seeing 


fa 3.0" 


2.0" 


« 2.5" 


Si 1.6" 


Candidates 


378 


329 


343 


391 


Targeted 


217 


45 


225 


219 


ELG redshifts 


44 


10 


84 


142 



emission lines above 8250A. These limit the identification of II/3 
and [OIII] above 2 « 0.65 — 0.7, however they do not interrupt 
identification of [Oil] and so the impact of the sky lines is limited. 

2.3 Galaxy redshifts 

We us e the 2dFGRS soft ware AUTOZ jCroom etal]|200ll) and 
RUNZ (ICollessetalJl200lh to search for emission features with 
which to identify galaxies in our observed sample and determine 
redshifts. AUTOZ performs an initial identification of each spec- 
trum, fitting absorption and emission features. Each fibre spectrum 
was then evaluated by eye to assign a redshift and quality rating, 
Qop (which ranges from to 5 depending on the confidence of the 
identification). Only objects with qop ^ 3 were accepted as positive 
identifications. 

Examples of the spectra obtained with the AAOmega instru- 
ment are provided in figure [3] The spectra are all binned to a bin- 
width of ~ lOA and key emission and absorption features are 
marked. We also show the unbinned data for the [Oil] feature (in- 
sets), where it is evident that the doublet nature of the feature is 
marginally detectable at the observed resolution. Although data 
were obtained on both the blue and red arms, only spectra from the 
red arm are plotted here as there are few features useful for iden- 
tification in the blue wavelength range given the signal to noise of 
our data. The key emission features that facilitate the identification 
of these galaxies with short exposure times, i.e. [Oil], H/3 and the 
[OIII] doublet are all evident in these spectra. 

There were no positive identifications of galaxies with just ab- 
sorption line features and no emission lines in the spectroscopic 
observations. It is likely however that any absorption line galaxies 
targeted remained unidentified given the relatively short exposure 
times, which make it difficult to confidently identify absorption line 
features in the spectra. 

A summary of the numbers of emission line galaxies identi- 
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Figure 3. Example spectra taken on the AAOmega spectrograph with the 385R grism, binned to lOA bins. Wavelengths of galaxy emission and absorption 
features are marked, however the features of key use in identification were the [Oil], H/3 and [OIII] doublet emission lines. The insets each show the [Oil] 
feature in close-up and unbinned, showing the doublet nature of the [Oil] feature to be marginally discernible given the resolution of the spectrograph. The 
red dashed-lines show the expected positions of the doublet peaks at rest-wavelengths of 3726A and 3729A. Objects 1, 3, 5 and 6 are all marked in figure|6]as 
showing discrepancy between spectroscopic and photometric redshifts. 



fied in our four fields is provided in table [2] Our most successful 
field was the COMBO-17 SI 1 field in which we were able to target 
219 ELG candidates in seeing conditions of « 1.6". In this field we 
identified 121 of the 219 candidates as being emission line galaxies 
from their emission lines with a confidence of qop ^ 3. In total we 
were able to identify 311 emission line galaxies over a combined 
area of 12.4 deg^, giving an average sky density of 25 deg^^. How- 
ever, three of the observed fields suffered poor seeing conditions 
of 2.0", limiting our ability to successfully identify objects in 
these fields. At worst, completeness was reduced to < 25% in the 
COSMOS field due to the seeing of ~ 3.0". However, in the more 
reasonable observing conditions encountered with the observations 
in the SI 1 field (where the seeing was 1.6") we find that the iden- 



tification rate is a more promising ~ 65%, with a sky density of 

^ AOdeg''^. 

Figure |4] gives the redshift distribution of the spectroscopi- 
cally confirmed galaxies in the Sll field. The plot incorporates 
all galaxies identified in the Sll field and the original photomet- 
ric redshift distribution from COMBO-17 data (black dashed line) 
also from the Sll field. Our spectroscopic sample follows the ex- 
pected distribution closely, with z^p^c = 0.66 ± 0.23 (compared 
to Zphot = 0.65 ± 0.21). There is some contamination from lower 
redshift (i.e. z < 0.5) galaxies and in the spectroscopic sample this 
is at a level of ~ 18% (compared to a level of ~ 23% obtained 
with the COMBO-17 photometric redshift sample). 

Figure |5] shows identifications as a function of source magni- 
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Figure 4. The redshift distribution of all successfully identified emission 
line galaxies from the four fields targeted using AAOmega is shown (red 
histogram). A total of 280 galaxies were successfully identified. The origi- 
nal photo-z redshift distribution incorporating all the galaxies selected using 
our z fa 0.7 colour cuts in the COMBO-17 Sll field is also shown (black 
dashed line). Both distributions are normalized to give an area under the 
histogram equal to unity. 



1 .0 
0.5 




20.6 



Figure 5. Number counts of objects observed as a function of SDSS i-band 
magnitude. The dark histogram shows counts of objects identified with 
emission lines (Af^m), whilst the pale histogram shows the total number 
of objects (A^r) observed in each magnitude bin. The top panel shows the 
fraction N^rn/^T ^s a function of i-band magnitude. Data are only shown 
from the Sll field as all other fields were limited by adverse seeing. 



tude in the Sll field. The lower panel shows number counts of spec- 
troscopically confirmed galaxies exhibiting emission lines (Nem, 
dark bars) and of all objects targeted with AAOmega fibres (Nt, 
pale bars), whilst the upper panel shows the fraction, Nf-TnlNr- 
The consistency of the 65% identification rate across our magni- 
tude range is evident and we are clearly reaching the i — 20.5 mag- 
nitude limit successfully. A small fall-off in the fraction of ELGs 
identified is observed in the fainter magnitude bins, however num- 
bers still remain high. 



0.2 0.4 

AAOmega emission 



0.6 0.8 
ine redshift 



Figure 6. Comparison of the spectroscopic redshifts of 24 galaxies in the 
central 0.5° X 0.5° of the Sll field with photometric redshifts from the 
COMBO-17 survey. The error bars show the az ~ 0.03 eiTor quoted by 
the COMBO-17 team for their photometric redshifts. The four numbered 
points (1, 3, 5 and 6) are objects in which the photometric and spectroscopic 
redshifts disagree by more than Sctz (the numbers refer to the spectrum 
numbers from figure[3). 



In figure (6] we compare our spectroscopically determined 
redshifts against the COMBO-17 photometric redshifts for those 
galaxies lying in the central 0.5° x 0.5° region covered by the 
COMBO-17 data. The vertical error bars represent the ct^ ~ 0.03 
error quoted by the COMBO-17 team for the photometric red- 
shifts. We find a total of 24 objects that have both COMBO-17 
photometric redshifts and spectroscopic redshifts from this work. 
Overall there appears to be good agreement between the data with 
just 4 outliers (taken here as a difference between the photomet- 
ric and spectroscopic results of Str^) having significantly different 
redshifts. The spectra for all four of these objects are given in fig- 
ure[3]and each of the outliers are marked in figure|6]by the spectrum 
number (1, 3, 5 and 6) from figure|3] We find a mean offset between 
the spectroscopic and photometric redshifts of Az = 0.01 ± 0.04 
(after excluding points 1 and 5). 

We show in figure|7]the distribution of spectroscopically con- 
firmed z > 0.5 ELGs (filled blue cicles), z < 0.5 ELGs (green 
triangles) and unidentified objects (red crosses) in the g-r vs r-i 
colour plane. It is evident that the z > 0.5 emission line galax- 
ies are reasonably evenly spread in the g-r vs r-i colour plane as 
are the objects without any discernible emission, although there is 
some bias in these to be towards the redder end of the selection in 
both r — i and g — r. The z < 0.5 emission line galaxies appear 
to be biased towards the upper left limits of the selection region, 
towards the low-redshift main sequence. These may be further re- 
duced by altering our constraints, however this would also remove a 
significant number of 2 > 0.5 objects at the same time. The model 
evolution track from figure[T]is again plotted for reference. 

Now looking at the properties of the galaxy spectra, we mea- 
sure the equivalent widths of the nebular emission lines by fitting 
Gaussian curves to the [Oil] 3727A, H/3 and [OIII] 5007A lines. 
We were able to measure equivalent widths with confidence for 
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Figure 7. Spectroscopic results from the Sll and e04 fields. We show 
objects identified as 2 ^ 0.5 emission line galaxies (filled blue circles), 
z ^ 0.5 emission line galaxies (green triangles) and objects with no identi- 
fied redshift (red crosses). The same evolution track as plotted in figure[T]is 
also shown. 



[Oil] 3727A, H/3 and [OIII] 5007A in 109, 53 and 51 of the galax- 
ies in our sample respectively. From these we determined mean 
equivalent widths of 23.0A, 8.12A and 8.98A for [Oil] 3727 A, H/3 
and [OIII] 5OO7A respectively. These mean equivalent widths are 
broadly consistent wit h other measurements of emissio n lines in 
late-type galaxies (e.g. lKennicut^ll992l : [shi et alJl2006h . In 27 of 
these galaxies we were able to measure all three of the above neb- 
ular emission lines with confidence and have attempted to evalu- 
ate the presence of AGN in our sample using the "blue diagnos- 
tic" constraints of Lama reille et al.l j2004h . which are based on the 
[OII]A3727/H/3 and [OIII]A5007/H/3 line ratios. This is shown 
in figure [S] where the solid line marks the estimated division be- 
tween AGN (above) and star-forming galaxies (below). The dashed 
lines mark the region of uncertainty between the two populations. 
In all, 22 of this sub-sample fall within the star-forming galaxy 
region, whilst the remaining 5 (2 of which have large uncertain- 
ties) fall within the uncertain region and none lie in the AGN re- 
gion. Within the reliability of the blue diagnostic diagram, we can 
say that our sample is dominated by star-forming galaxies and this 
method shows no positive evidence for any AGN contamination of 
our sample although there are a small number of borderline cases. 

Figure [9] shows a composite spectrum of all of the confirmed 
emission line galaxies over all redshifts, with significant emission 
and absorption features labelled. The key emission lines used in our 
spectral identification (i.e. [Oil], H/3 and [OIII]) are clearly evident. 
We also see the Balmer absorption features red-wards of the [Oil] 
emission, whilst the weak ELG 4000A break is also apparent in this 
composite. 



3 CLUSTERING 

3.1 Angular Correlation Function 

We now evaluate the angular correlation function for a sample of 
galaxies selected based on our three photometric selections. The 
datasets taken from SDSS DR6, as described in the previous sec- 




[OII]Ai727/H^ (EW) 



Figure 8. "Blue diagnostic" diagram based on ^Lamareill e et al.l j2004h . 
Line ratios are plotted for the sub-sample of our spectroscopically observed 
sample for which we have equivalent widths for the [OIIJA3727, H/3 and 
[OIIIJA5 007 nebular emission li nes. The solid line marks the limit esti- 
mated bv lLamareille et alj ilOOj) between star-forming galaxies and AGN 
and the dashed lines show the region of uncertainty. In total 22 objects lie 
within the star-forming region of the diagnostic plot and are marked by 
filled blue curies. A further 5 lie within the overlap region (cyan stars) and 
none of the objects lie within the AGN region. 



tion, are used for this purpose. We calculate the angular-correlation 
function of the samples using the estimator; 



w{e) 



ur dd 

rw^R 



(1) 



where DD and nn are the numbers of galaxy-galaxy pairs and the 
total number of galaxies respectively. For these calculations we use 
a random catalogue which exactly matches the sky coverage of our 
SDSS galaxy samples and with a factor of 20 more random points 
than galaxies in each of our galaxy samples. The total number of 
random points is given by ur and DR is simply the number of 
galaxy-random pairs. Statistical errors are estimated using field-to- 
field errors, using 16 separate fields within our complete field. Our 
results for the three photometric samples are shown in figure [Tol 
where the blue triangles, green squares and red crosses show the 
low, mid and high redshift samples respectively. For comparison 
we also plot w{9) for a photometrically selected L RG sample from 
SDSS data at z « 0.7 jSawangwit & Shanksl26o 9). The clustering 
amplitude of all three of our samples is significantly lower than the 
LRG amplitude. 

From these measurements of the angular correlation function, 
we now estimate the 2-point correlation function (2PCF, ^(r)) us- 
ing Limber's formula. We make an estimate of ^(r) for each of 
the samples using a double power-law with a central break, i.e. 
^(< rb) = (t^)-'" and ?(> n) = i^)-'' ■ This is then 
combined with our best estimate of the redshift distribution (based 
on the COMBO-17 photometric redshift data for the low and mid 
redshift samples and the spectroscopic redshift data for the high 
redshift sample) to calculate the resultant w{8) wit h Limber's for- 
mula. A full treatment of this calculation is given bv lPhillipps et al.l 
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Figure 9. Composite spectium of the 280 successfully identified emission line galaxies. The key emission line features used for identification are clearly 
visible: [Oil], H/3 and [OIII], whilst absorption features which are difficult to observe in individual spectra are now evident. 



( Il978h . We then perform a fitting, in the range 2' < 6 < 20' 
to our data. The best fitting models are plotted with the data in fig- 
ure [Tol whilst the associated parameters are listed in table [3] We 
find reasonable fits to both the low- and mid-redshift samples, the 
low-redshift sample being well fitted by a double power law with a 
break at 0.5/i~^Mpc and the mid-redshift sample by just a single 
power-law. We note however, that we struggle to fit to the high- 
redshift sample with either a double or single power law. This is 
largely due to strong deviations from a simple power law trend at 
separations of < 2'. This results in large values for our attempts 
to fit the correlation function at in this range. The angular correla- 
tion function does however return to a simple power law at separa- 
tions of 2' < ^ < 20' where we are able to provide a reasonable 
power-law fit using the Limber method. 

From our estimates of ^(r), we now go on to estimate the bias 
of each sample. The biasing parameter, b, quantifies the relative 
clustering of a given galaxy popula tion compared to th e underly- 
ing dark matter distribution ( Tegmark & Peeble3ll998[) . This can 
be expressed as the following: 

igai{r) = b'^DMir) (2) 

Here, ^gai{r) is the 2PCF of the galaxy sample and (,DM{r) 
is the 2PCF of dark matter at the same epoch. We determine the 
dark matter correlation function by first using the CAMB (Code 
for Anisotropics in the Microwave Background. iLewis et al. 20001 
software to estimate the DM power spectrum at the mean redshifts 
of each of our galaxy sample s. The power sp ectrum is calculated 
using the HALOFIT model dSmith et al.l2003l) to fit non-linear fea- 
tures, at each of the mean redshifts of our samples. With the DM 
power spectra calculated at each redshift we then simply calculate 
the corresponding 2PCFs via the Fourier transform: 

CdmW = 4n r k^P{k)^^^^dk (3) 
Jo kr 

We now estimate the bias by evaluating the dar k matter and 
the galaxy 2PCFs to a maximum s eparation of 20Mpc jCroom et al.l 
l2005l : ISawangwit & Shanksll2009l) . This limit, restricts the calcula- 
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Figure 10. The angular correlation functions, w{6) for our three photo- 
metric redshift selections. Blue crosses, green diamonds and red squares 
represent the low- (z < 0.4), mid- (0.4 < z < 0.6) and high-redshift 
(z > 0.6) samples respectively. The best fitting power law models are 
plotted through each s et of d ata. The dashed black line shows the w{8) of 
ISawangwit & Shanks! J2009l) for their sample of photometrically selected 
Luminous Red Galaxies at 2 = 0.68 ± 0.1. 

tions to the linear regime at which our fits to the correlation function 
are still valid. Thus the biasing parameter can be estimated using: 

Idm(20) 

where ^{x) is given by: 

ax) = r r\{r)dr (5) 
Jo 

We show ^(20) for each of our three redshift samples in fig- 
ure [TT] (denoted by the stars). For comparison we also plot ^(20) 
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Table 3. Comoving correlation lengths, ro and power-law slopes, 7, for the double power-law model used to provide fits to the angular correlation functions 
for each redshift selection. 



z 


T-f, (h ^Mpc) 


'■o(< n) (ft-iMpc) 


7(< rb) 


'■o(> rt) (h-iMpc) 


7(> rb) 


b 


0.29 ±0.05 
0.44 ±0.08 
0.65 ±0.21 


0.5 
n/a 
0.5 


'^■•^"-0.03 

n/a 
r, on+0 05 


r, 91 +0.02 

^■^^-0.01 
n/a 

9 79+0.01 
^- '^-0.01 


9 f;r+0.05 

q f.9+0.06 
'^■°'^-0.06 
c- 00+0.12 
b.8«_0.i2 


1 K4+O.02 
^■3^-0.02 
1 65+°°l 
1 OO+0.03 
^■»J_0.01 


0.70 
0.92 
1.46 
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Figure 11. Estimated ^(20) plotted versus redshift for each of our three 
photometric samples (st ars). Also plotted is th e |(20) for late type galaxies 
from the 2dFGRS from lNorberg et alj j2QQ2h , split by absolute magnitude 
(triangles). The dashed lines show the long-lived model normalised to the 
ELG data-points, whilst the dash-dot lines show the stable clustering model 
also normalized to the ELG data. The solid lines project the clustering of 
each population with no evolution. 



for the 2dFGRS late type galaxy samples o f iNorberg et all j2002h . 
These are split into absolute magnitude bins of —18 < Mbj — 
5logio{h) < -19, -19 < Mbj - Uogw(h) < -20, -20 < 
Mbj-5logio{h) < -21 and -20.5 < Ahj -5logio{h) < -21.5 
and are calculated based on the correlation parameters given in their 
table[3] Based on the redshift and apparent magnitude distributions 
of our three samples, we estimate absolute magnitude ranges of 
Mb^ - blogioih) = -19.9 ± 0.3, -20.0 ± 0.3 and -20.9 ± 0.4 
for the low, mid and high redshift samples respectively. These esti- 
mates include K+e corrections based on the r = 9Gyr SFR model 
given in section |2] 

Also plotted in figure [TT] are three simple clustering evolution 
models: the long-lived model (dashed lines), stable clustering (dot- 
dashed lines) and no evolution of the comoving-space clustering 
(solid lines). All three models have been normalised to each of the 
ELG clustering amplitudes. 

The long-lived model is equivalent to assuming that the galax- 
ies have ages of order the Hubble time. The clustering evolution 
is then governed by their mot ion within the gravitational potential 
( lF^ll996l : lte>om et al.ll2005l) . The bias evolution is thus governed 
by: 



b{z) 



1 + 



m - 1 

Diz) 



(6) 



where D{z) is the linear growth rate and is determined using the 
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fitting formulae o f lCarroll et all ( 1 19921) . For the bias, b, we use the 
values given in tabled 0.70, 0.92 and 1.46 for the low, mid and 
high redshift samples respectively. The stable clustering model rep- 
resents the evol ution of virialised structures and is characterised by 
jPeacocldl 19991) : 



OC (1 + 2 



,7-3 



(V) 



where r is the comoving distance and 7 is the slope of the clus- 
tering correlation function. Finally, the no evolution model simply 
assumes that there is no evolution of the clustering in comoving 
coordinates. 

Before analysing this data, it is important to note that the var- 
ious samples have not been refined to be precisely equivalent in 
that neither the number densities or absolute magnitudes have been 
ideally matched between any of the samples. Further to this, the 
clustering strengths of the z ~ 0.3 and z ~ 0.45 are based on the 
photometric redshift distributions from the COMBO-17 Sll field, 
which as yet have not been fully calibrated. This comparison is 
intended to give only an indication of the potential clustering evo- 
lution of the the samples discussed. 

Looking at the low and mid-redshift ranges we find that 
the clustering measurements suggest an increase in the cluster- 
ing amplitudes of the emission line galaxies from redshifts of 
2 = 0.3 — 0.5 to more recent epochs if they are to evolve into 
the 2dF samples with comparable absolute magnitude ranges (i.e. 
the -19.5 > Mb > -20.5 2dF sample). 

The high-redshift sample appears more consistent with a no 
evolution model if it is to evolve into the most comparable 2dF 
samples (i.e. —20.5 > Mb > —21.5). However, the stable and 
long-lived clustering evolution models are still only 1 — 1.5(t above 
the brightest 2dF ^ (20) value and therefore cannot be ruled out by 
the present accuracy of the 2dF late-type data. 



3.2 Redshift-space Correlation Function 

We now estimate the redshift-space correlation function, (,{s), us- 
ing 2 ^ 0.5 galaxies identified with qop ^ 3 from the four fields 
observed with AAOmega. The redshift distribution from figure |4] 
was used to create random catalogues with which to perform the 
auto-correlation analysis. In each field we use a catalogue of 20 x 
the number of random points as galaxies in that field. In total this 
calculation encompasses 276 galaxies across 12.6deg^. 

We use the correlation estimator given in equation [T] to deter- 
mine 5(s), whilst errors are estimated using Poisson errors. The 
result is shown in figure [T2I (filled square points). We fit the ^(s) 
measurement with a single power-law (noting that the break used in 
the double-power laws previously lies below the range of our ^(s) 
result) and find a best-fit (using a fixed slope of 7 = 1.8) given by 
a clustering length of So = 6.4 ± 0.8h~^ AIpc (solid line). 

We also show the expected ^(s) determined from the power- 
law form of ^(r) given by our estimate of w{6) (dashed line). To 
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Figure 12. The redshift-space correlation function (^(s)) for tlie full sam- 
ple of spectroscopically identified objects. The data-points show the cor- 
relation function of the spectroscopic sample of galaxies incoiporating the 
Sll, COSMOS, e04 and d05 fields. The solid line shows a best fit to the 
data points (with a fixed slope of 7 = 1.8), which is characterised by a 
correlation length of So = 6.4 it The dashed line is the ^(s) 

determined from the angular correlation function estimate of the real-space 
congelation function (ro = 5.88/i~^ Afpc, 7 = 1.83 at r > lh~^ Mpc). 
This incorporates the coherent infall (/3 = 0.54) and random pairwise ve- 
locity (a = 500A:ms~^) effects. 



CMB. Further to this, the large scale nature of the effect lends it- 
self well to the use of photometrically selected galaxy populations 
as detailed r edshift information is unnecessary. In this vein, several 
authors (e.g.lFosalba et alj|2003l : IScranton et al]|2003l : ICabre et al.1 
l2006l : lRassatetal.ll2007h have used a number of galaxy samples to 
attempt measurements of the ISW signal in the WMAP 1st and 3rd 
year data. The samples used have mostly been photometrically se- 
lected LRGs at redshifts of z < 0.6 and simple magnitude selected 
samples. 

We now use our ELG sample to attempt to measure the ISW 
effect in WMAP 5-year data at our sample redshifts of 2 ~ 0.3, 
z ~ 0.5 and z ~ 0.7. As stated much of the ISW work done thus 
far has been with LRGs and magnitude cut samples at z < 0.6. 
Our use of the ELG samples provides the benefit of extending to 
greater redshifts, whilst also using an alternative galaxy popula- 
tion. This in itself has benefits and draw-backs. Firstly, the mea- 
sured signal will be heavily dependant on the bias of the sample 
(i.e. how well the sample traces the dark matter structure and hence 
the gravitational potential). Given that the ELGs are less clustered 
than the LRG samples used thus far, we therefore expect to mea- 
sure a weaker signal, making the measurement potentially more 
difficult. The potential gain in the low clustering strengths of the 
ELG samples however, is that they are less likely to reside in rich 
clusters and so we may expect the ISW signal to be less affected 
by the SZ effect produced as CMB photons pass through hot intra- 
cluster gas. This potentially provides an interesting alternative to 
the highly clustered LRG samples used in a number of the previous 
studies. 



do this, we take the power-law form and apply both coherent infall 
and random pairwise velocity effects. The coherent infall imprint 
on the correlation function is characterised by the infall parameter, 
/3, which is given by; 

n (z)°-^'^ 
P = "'Y (8) 


where Q.m{z) is the mas s density at t he required redshift and b 
is again the sample bias ( lKaiseilll987l) . Using the value of 6 = 
1.46 from table [3] this gives a value of the infall parameter of 
/3 = 0.54. For the random pairwise velocities we use a value of 
a = 500fcms~^ as a reason able estimate of the r andom motions 
based on the 2dFGRS data faawkins et alj|2003l) . The expected 
f (s) w as then calculated using the relations given in'Haw kins et al.l 
Based on this calculation, we see from figure[T2]that we find 
good agreement between the ^(s) of our observed spectroscopic 
sample and that determined from the photometric data, within the 
associated errors. 



4.2 Data 

For this cross-coiTelation, we have used the W, V, Q and K 
band temperature m aps from the WMAP 5-year data release 
jHinshaw et al.ll2008h . We use the full-resolution maps in all cases. 
Before performing the cross-correlation we apply two masks to the 
data. The first is the WMAP KPO mask (Bennett et al. 2003) which 
removes the majority of the galactic (Milky Way) foreground and 
is the most rigorous mask provided by the WMAP team. Secondly 
we mask the data to match the coverage of our SDSS DR6 galaxy 
samples, which is described further below. 

Pixelised sky-density maps are constructed from each of the 
three galaxy samples using the HEALPIX software. These are con- 
structed with a resolution identical to the WMAP temperature maps 
characterised by the HEALPIX parameter NSIDE = 512 (pixel 
width~ 7'). We then limit our galaxy sample to incorporate only 
the contiguous north galactic pole region of the SDSS. Thus our 
sample is limited to 100° < R.A. < 270° and stripes 39, 42 and 
43 are also excluded. 



4 INTEGRATED SACHS WOLFE EFFECT 
4.1 Overview 

As described in the introduction the ISW effect is characterised by 
the energy boost that CMB photons experience as they cross tempo- 
rally evolving gravitational potential wells in an accelerating Uni- 
verse. The effect is therefore a potential tool for placing constraints 
on the acceleration of the Universe as characterised by the cosmo- 
logical constant, A. Given the effect's large scale, the production of 
all-sky CMB maps from the WMAP experiment has made it possi- 
ble to attempt measures of the effect through the cross-correlation 
of galaxies (as tracers of the large gravitational potentials) with the 



4.3 Method 

Following the work of iFosalba et"d] ( l2003l) and lCabre et"ZI ( l2006h , 
we use the estimate the cross-correlation of the galaxy and WMAP 
data as the expectation value of the product of the galaxy over- 



density, Sg 



and the normalised CMB anisotropy temper- 



ature, AT = T — T as a function of the angular separation, 6. This 
is given by: 



WTg{e) = 



(9) 
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Again following lFosalba et alj ( l2003h and lCabre etal] j2006l) . 
the form of the ISW as probed by a given galaxy population can be 
expressed by the following Legendre polynomial expansion: 



^r'J^ie) - E ^P<(cose)cFT^(0 (10) 
I 

(l) is simply the ISW/galaxy population power spec- 
trum as given by: 



H{2 



(2/ + 1)2 7 Wisw{z)WG{z)^^P{k)dz (11) 

where P(k) is the mass power spectrum and Wisw{z) and Wg{z) 
are given by: 



Wiswiz) = 3n„ 



HoY d[D{z)/a] 
dz 



(12) 




100 
6 (arcmins) 

Figure 13. Cross correlation between the low-redshift galaxy sample and 
the WMAP V band data. The solid line shows the predicted model. Errors 
are field-to-field based on splitting the data sample into 16 distinct seg- 
ments. 



Wg{z) = b{z)(t>{z)D{z) 



(13) 



where D{z) is the linear growth rate and b{z) is the bias of the 
galaxy population (taken from section [5}. 4>{z) is the galaxy selec- 
tion function, set from the n{z) distribution of each of the galaxy 
samples. 



4.4 Results & Error Analysis 

We perform the cr oss correlation u sing the NPT (N-point spatial 
statistic) software i Grav et al]|2004h with the weighting for each 
pixel given by the galaxy density, Sg and the CMB anisotropy tem- 
perature, AT. The results are shown in figures [Ts] to [15] for four 
WMAP bands: W, V, Q and K. We als o plot the predicted result 
using predictions based on equation 6 of lCabre et al.l j2006l) . 

We estimate the errors on the cross-correlation analysis using 
field-to-field errors. For this purpose, we split the studied region 
into 16 approximately equal area sub-fields and re-calculate the 
cross-correlation within each sub-field. The error in each angular 
bin is thus estimated as a j ^pn, where a is the standard deviation 
across the sub-fields and n is the number of sub-fields. 

Summing over all bins at ^ < 100' we find amplitudes 
for mTg(< 100') in the WMAP W-band of 0.25 ± 0.27fiK, 
0. 17 ± 0.20^4 A' and 0.17± O.lB/iii' for the low, mid and high red- 
shift samples respectively. Similar results are obtained with the V 
and Q bands, whilst the K-band (which has a greater level of galac- 
tic contamination and a lower resolution) is less consistent, giving 
signals ofuiTgiK 100') = 0.13 ± 0.36^*^:, -0.16 ± 0.29^iC and 
0.38 ± OASfiK for the low, mid and high redshift samples respec- 
tively. 

We also evaluate the significance of the observed correlation 
by repeating the cross-correlation with rotated realisations of the 
WMAP data. This method uses the data itself in place of random 
realisations by rotating the masked WMAP data in 30° steps in 
galactic longitude. We note at this point that rotating in R.A. would 
lead to the galactic plane entering the field of view and although 
the galactic plane region is masked, it would reduce the number 
pixels in the analysis significantly. For consistency, we also rotate 
the WMAP KpO mask before applying it to the galaxy density map. 
The result of this treatment, using the high redshift sample, is given 
in figure [T6l Here the iUTg(< 100') signal is plotted as a function 
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Figure 14. As in figure [T3] but with our mid-redshift sample of emission 
line galaxies. 



of rotation of the WMAP data through a full 360° in galactic lon- 
gitude. The dotted line shows the non-rotated signal. Again we see 
that the positive signal that we see in the data does not appear statis- 
tically significant, with the rotated results showing a large amount 
of scatter around WTg{< 100') = O^K and two of the results (at 
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Figure 15. As in figure [T3] but with our high-redshift sample of emission 
line galaxies. 
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Figure 16. Cross-correlation signal, WTg[< 100'), for the high redshift 
galaxy sample as a function of rotation of the WMAP 5yr data in galac- 
tic longitude. The plotted errors are field-to-field en'ors based on the seg- 
mentation of the data into 1 6 distinct regions and the dotted line shows the 
measurement from the non-rotated data. 
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Figure 18. Cross-correlation signal, WTg[< 100'), for the combined low-, 
mid- and high-redshift galaxy sample as a function of rotation of the WMAP 
Syr data in galactic longitude. The plotted eiTors are field-to-field errors 
based on the segmentation of the data into 3 X 16 distinct regions. Again 
the result from the non-rotated data is shown by the dotted hne. 
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Figure 17. Cross correlation result averaged across the 3 X 16 red- 
shift/segment samples. 

30° and 90°) showing more significant positive correlation than the 
non-rotated result. 

We now attempt to improve our statistics by combining the 
low, mid and high redshift results. Figure[T7]shows the mean of the 
3 X 16 separate cross-correlation results. The errors are again given 
by the field-to-field errors, this time across the whole 48 sample set. 
Again we see a positive signal that appears to show some agreement 
with the model. In the W, V, Q and K WMAP bands we derive 
signals of uiT9(< 100') = {Q.2Q ±G.12)^K, (0.20 ± 0.12)^A', 
(0.18 ± Q.12)^J.K and (0.11 ± Q.16)^J.K respectively. Repeating 
the rotation analysis (figurefTSt. but with the combined sample, we 
again find significant scatter about WTg{< 100') — 0. Indeed, a 
stronger signal is again found at some rotation angles than at the 

zero position. 

Comparing this analysis to previous results, ICabre et^ 

( I2OO61) obtained a signal at 100' of WTg ~ 0.6 ± 0.3iJ,K using 
a sample of 2 ~ 0.5 LRGs and a signal of wtq =~ 0.65 ± 0.2fiK 
with a 20 < r < 21 magnitude selected sample with a median 
redshift of 5 = 0.28. In addition, Fosalba & Gaztaiiag^ j2004l) 
claimed a detection of 0.35 ± O.lAfiK at angular scales of 6 = 



4° — 10° using APM galaxies, but found their signal was dom- 
inated by the SZ component at scales of < 4°. Our observed 
signal is consistent with the weaker signal expected of the ELG 
samples, based on both the model predictions and comparison with 
correlations based on more strongly clustered populations. How- 
ever, the estimated statistical and particularly the systematic errors 
mean that no detection of the ISW effect can be claimed at this 
point. 



5 CONCLUSIONS 

Colour selected samples are an extremely useful tool in modern as- 
tronomy and cosmology. They offer a cheap route to large galaxy 
redshift surveys, facilitating investigations in several key areas of 
interest, including the study of Dark Energy via both the Baryon 
Acoustic Oscillations and the Integrated Sachs Wolfe. With this 
in mind we have developed photometric selections with which to 
identify galaxies in three broad redshift ranges characterised by 
z = 0.29 ± 0.05, 5 = 0.44 ± 0.08 and i = 0.65 ±0.21. Applying 
these to SDSS data, we are able to select galaxy samples covering 
the given redshift ranges down to magnitudes of i ~ 20.5. 

The calibration of the high-redshift sample has been per- 
formed, using the AAOmega spectrograph at the AAT to provide 
spectroscopic redshifts for this sample. From these observations, 
we have shown that it is possible to target and spectroscopically 
identify z > 0.5 star-forming galaxies with integration times of 
just 1 hour on the 4m AAT. The results of this have shown our 
high-redshift selection to work well, giving a close match to the 
expected redshift distribution, with z — 0.66 ± 0.23. Further, this 
has proven that our selection can be a firm basis on which to con- 
duct large scale spectroscopic surveys of 2 > 0.5 emission line 
galaxies. 

We have investigated the clustering properties of our three 
photometric samples using the angular correlation function. By fit- 
ting the angular correlation function using Limber's formula, we 
estimate the real-space clustering properties of the samples and 
find that both the low and high redshift populations are best fit 
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by double power laws, whilst the mid-redshift population seems 
best fit by a single power-law. From these calculations, we estimate 
clustering lengths (at r > 0.5h'^Mpc) of 2.65 ± O.Uh'^Mpc, 
3.62 ± 0.12h-^Mpc and 5.8 ± 0.12h-^Mpc for the low, mid 
and high redshift samples respectively. Further to this, from our 
spectroscopic observations we measure a clustering length for the 
high redshift sample of 6.4 ± 0.8h~^ Mpc in agreement with the 
measurement based on the angular clustering measurement. Com- 
paring these clustering measurements with comparable populations 
of late-type galaxies from 2dF data, we note that the high redshift 
sample appears to have an unexpectedly high clustering strength. 

The development of photometric redshift selections has a 
number of scientific applications, not least the evaluation of the 
ISW effect. We have used our photometric selection to evaluate 
the ISW effect in the region of the SDSS by cross correlating 
the density fluctuations in our galaxy distributions with the CMB 
anisotropies from the WMAP 5yr data. The results obtained using 
all three datasets show a positive correlation in accordance with 
the predicted model. However, none of these three prove to be sig- 
nificant, with signals in the WMAP W-band of (0.25 ± 0.27)^iK, 
(0.17 ± 0.20)^iA' and (0.17 ± 0.16)fiK for the low, mid and high 
redshift samples respectively. We attempt to improve the statistics 
by combining the three redshift samples, which results in a signal 
of (0.20 ± 0.12) ^iK when cross-correlated with the WMAP W- 
band, still only marginal (1.67cr). Also in tests of contamination 
by systematics, we found similar results at arbitrary angles of ro- 
tation between the CMB data and the ELG samples which means 
that no detection of the ISW can be claimed above the random and 
systematic noise in this analysis. 
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